The fast spiral-in of the companion star 
to the black hole XTE J1118+480 1 
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ABSTRACT 

We report the detection of an orbital period decay of P = — 1.83 ± 0.66 ms yr _1 in the black 
hole X-ray binary XTE J1118+480. This corresponds to a period change of —0.85 ± 0.30 [is 
per orbital cycle, which is ~ 150 times larger than expected from the emission of gravitational 
waves. These observations cannot be reproduced by conventional models of magnetic braking 
even when including significant mass loss from the system. The spiral-in of the star is either 
driven by magnetic braking under extremely high magnetic fields in the secondary star or by a 
currently unknown process, which will have an impact on the evolution and lifetime of black hole 
X-ray binaries. 



Subject headings: black holes: physics — gravitation 
magnetic field — X-rays: binaries 

1. Introduction 

Angular momentum loss (AML) in short- 
period black hole X-ray bin aries are driven essen - 
tially by magnetic braking (IVerbunt et al.lll98ll). 
gravitational radiation ()Landau fe Lifshit zlTl 962; 



stars: individual XTE J1118+480 



stars: 



Tavlor fc Weisberdll982l) . and mass loss (jRappaport et al 



19821 ) , with possibl e addi tional contributions from 
jets (jKing fe Kola Il999f ) . Measurements of or- 



bital period variations with time constrain the 
role of t hese processes in the evolution of such 
binaries (|Verbuntl ll993). 

Magnetic braking is proposed as the main 

mech anism for AML in compact binaries (jRappaport et al 
1983). It produces the shrinking of the binary or- 
bit and maintains the donor star in contact with 
its Roche lobe, therefore, sustaining stable mass 
transfer. Although this is an essential mechanism 



to understand compact binary evolution its ac- 
tual prescri ption is current ly debated and not well 
established ( Ivanova l2006h . 

XTE J1118+480 is one of only 18 galactic X- 
ray binaries which contain a dynamical black 
hole ( Casaresl 2007 ). With an orbital period of 
4.1 hours it is the most compact black hole X-ray 
binary known. Determination of the times at the 
inferior conjunction of the ~ 0.2 M Q secondary 
star with respect to the ~ 8 M Q black hole have 
been obtained at different epochs since its discov- 
ery on UT 20 00 March 29 by the Ro ssi X-ray Tim- 
ing Explorer (jRemillard et al ] l2000l) . Previous at- 
tempts to estimate the orbital period derivative 
failed due to the sma ll baseline of the observa- 
tions ( Johannsen 20091 ). 
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2. Observations 

We have conducted new spectroscopic observa- 
tions of XTE J1118+480 using the 10.4m Gran 
Telescopio Canarias ( GTC) equipped with the 
OSIRIS spectrograph (|Cepa et al.ll2000l 120031) at 
the Observatorio del Roque de los Muchachos 
in La Palma (Canary Islands, Spain). Ninety 
seven medium-resolution spectra (X/SX ~ 2,500) 
were obtained on UT 2011 January 7, Febru- 
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ary 8, and April 25-25, 36 and 36 spectra in 
each night respectively-. In Fig. [1] we display 
the radial velocities of the secondary star in its 
orbital motion around the center-of-mass of the 
system, derived using the cross-correlation tech- 
nique with a template stellar spectrum prop- 
erly broadened with a rotational velocity of 
Dsini = 100 kms -1 ([Gonzalez Hernandez et al 



l2006l l2008h . The radial velocity points, which 
spread over ~ 3 and a half months, provide a 
new determination of the current orbital period 
of P orb = 0.16993379 ± 0.00000047 d, which is 
smaller although still consistent with the orbital 
period measureme nt previously determ ined on UT 
2000 December 1 (jTorres et al.ll2004) . In Tabled] 
we list these orbital period measurements, P OT b, 
the updated dynamical masses of the black hole, 
MbH) an d the secondary star, Af 2 , the mass ra- 
tio, q, the rotational velocity, usini, the orbital 
inclination, i, the orbital semiamplitude velocity, 
k>2, the mass function, f(M), the current radius 
of the secondary star, R2, and orbital separation 
of the system components, a c . 

3. Orbital period decay 

The spectroscopic data were used to derive 
three new times, T„, of the inferior conjunction 
of the secondary star in this system (see Table [2]) . 
Assuming a constant rate of change of the orbital 
period, the time, T n , of the nth orbital cycle can 
be expressed as T n = T + P n n+ ^P Pn 2 , where Pq 
is the orbital period at time To of the reference cy- 
cle (n = 0), P is the orbital period time derivative, 
and n, the orbital cycle number. We use the IDL 
routine CURVEFIT and obtain T = 2451868. 8921± 
0.0002 d, P = 0.16993404 ± 0.00000005 d, and a 
period derivative of P = -(5.8 ± 2.1) x 10~ u s/s 
with a reduced xl = 1-7 with v = 3. A linear 
fit (P = 0), and a third-order polynomial fit (in- 
cluding P), provide worse fits with \t = 2.9 and 
2.3, respectively. In Fig. |2]we have depicted the 
orbital phase shift, defined as <f> n = T "^ T " — n, 
of each of the T n values as a function of the or- 
bital cycle number n, together with the best-fit 
second-order solution. This figure shows a clear 
deviation from the null variation and that P is 
negative. Our result, which can be expressed as 
P = —1.83 ± 0.66 ms yr _1 , represents the first de- 
termination of the orbital shrinkage in a low-mass 
black hole X-ray binary. 
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Fig. 1. Top panel: radial velocities of the secondary star 

in the black hole X-ray binary XTE J1118+480 obtained from 
the GTC/OSIRIS spectroscopic data taken on the three nights 
of UT 2011 January 7 (green filled circles), UT 2011 Febru- 
ary 8 (blue filled diamonds) and UT 2011 April 25 (red filled 
triangles), folded on the best-fitting orbital solution. Bottom 
panel: residuals of the fit, with a rms of ~ 20km s _1 . 
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Fig. 2 . Top panel: orbital phase shift at the time of the 

inferior conjunction (orbital phase 0), T n . of the secondary star 
in the low-mass black hole X-ray binary XTE J111S+480 ver- 
sus the orbital cycle number, n, folded on the best-fit parabolic 
fit. Green filled circles are spectroscopic determinations, the 
blue filled triangle is a photometric measurement and red di- 
amonds are the new GTC/OSIRIS spectroscopic determina- 
tions. Solid error bars show the uncertainties associated to the 
T n determinations, whereas the dashed error bars also include 
the uncertainties associated to the To and Pq determinations. 
Bottom panel: residuals of the fit of the T n values versus the 
cycle number n. 
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Table 1 

klnematical and dynamical binary parameters of xte j1118+480. 



Parameter 


Value 


Reference 


v sin i 


lOOti! krns" 1 


[i] 


i 


68 ± 2 


[2] 




708.8 ± 1.4 km s -1 


[1] 


Poxb 


0.1699339(2) d 


m 


Porb 


0.1699338(5) d 


[4] 


q = M 2 /M BH 


0.027 ± 0.009 


[i] 


f(M) 


6.27 ± 0.04 M 


[i] 


M BH 


8.30±°;^ M 


[4,5] 


Mi 


0.22 ± 0.07 M 


[4,5] 


R 2 


0-37lo;°5 R 


[4,5] 


a c 


2-64+g °| R 


[4,5] 



Not e. — R eferences: [1 ] Gonzalez Hcrnande^^^i 1 
120061) ; [2] IGelino et ail 120061); 131 ITorres et a .1 
|2004j) ; [4] This work; [5] IGonzalez Hernandez et all 



Discussion and conclusions 



General Relativity (Mi sner et al.lll973 ) predicts 
a secular periastron precesion period of ~ 23 yr 
for XTE Jl 118+480. In the case of non-zero or- 
bital eccentricity, this would affect the period in- 
ferred from succesive inferior co njunctions (see 
e.g. eq. 20 in lPal fe Kocsisll2008l ). Using our data 
we set an upper-limit to the eccentricity of the 
black hole binary orbit of e < 0.0067 (95% confi- 
dence level), and conclude that relativistic perias- 
tron precesion could only explain up to 50% of the 
measured orbital period decay. 

The loss of energy through emission of quadrupole 
gravitational radiation, as predicted by general 
relativity, provides an orbital period derivative 
for this binary system of Pqr ~ — 0.01 ms yr^ 1 



(from eq. 8 in Taylor fc Weisberg] 1982) which is 



far below the observed rate of variation of the or- 
bital period. On the other hand, the presence of 
jets in XTE J1118+480 has also been postulated 
bas ed on X-ray and infrared observations (see 
e.g. iMiller et alj|2002j ). However, the consequen- 
tial AML due to j ets even for relativistic jets 
(|King fe Kolbl Il999h is < 1000 times lower than 
the AML due to magnetic braking. The presence 
of a circumbinary disk is also proposed as possi- 
ble AML mechanism in cataclysmic variables (see 
e.g. iTaam et al.ll2003l) . and t he detection of mid- 
IR 4. 5-8 /xm excess emission (jMuno fc Mauerhan 
20061 ) do support its existence in XTE Jl 118+480. 



However, the estimated total mass of this disk is 
too small to affect the orbital angular momen- 
tum of this X-ray binary. Therefore, these effects 
cannot account for any significant additional con- 
tribution to the apparently too large orbital period 
decay measured in this black hole binary. 

Most conventional magnetic braking and mass 
loss models for the binary parameters of XTE J1118+480 
predict a much lower orbital period decay than 
we measure (see Fig. [3]). Our observations ap- 
pear inconsistent with these models. Only a 
very restrictive, rather unplausible subset of mod- 
els, where nearly all the mass transferred by 
the secondary star is also lost by the system 
(f3 = —Mbh/M 2 ~ 0), may lead to values of the 
orbital period derivative close to our observations 
(see Fig. E]). 

Magnetic braking could be enhanced by anoma- 
lously high magnetic fields in the secondary star 



We ha ve estimated, from eq. 5 in lJustham et al 
(|2006l ) and assuming that the mass lost by wind is 
eq ual to the mas s trans fer rate (derived from eq. 9 



King fc Koibl (|l999T )V that a magnetic field at 
the surface of the companion star B s > 10 — 20 kG 
is needed to explain the observations. The angu- 
lar momentum loss due t o magnetic braking is de - 
termined from eq. 9 in Ijohannsen et al" ( 20091 ). 
This magnetic field is 1-2 orders of magnitude 
larger than typical magnetic fi eld strengths in 



highly-rotating low-mass stars (|Phan-Bao et al 



2009). However, we cannot discard that the donor 
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Table 2 

Time at inferior conjunction of the secondary star in XTE Jl 118+480. 



JV T n - 2450000 a ST n Reference 






1868.8916 


0.0004 


[1] 


0(5 


1880.1086 


0.0004 


[2] 


904 


2022.5122 b 


0.0004 


m 


6950 


3049.93347 


0.00007 


[4) 


21772 


5568.6936 


0.0003 


m 


21960 


5600.6413 


0.0002 


m 


22407 


5676.6017 


0.0002 


[5] 



a Timc at inferior conjunction, T^, of the sec- 
ondary star in the low-mass black hole X-ray bi- 
nary XTE J1118+480, and uncertainties, <5T„, cor- 
responding to a given orbital cycle number, n. 

b From photometric measurements. 



NOTE.— Refere n ces: [1] IWagner et al.l d2001j) : 
[21 ITorres et al.l ll2004|l: [3] "lZurita et al.l j2003li : 
[4] IGonzalez Hernandez et aU d2008D : [51 This work. 




Fig. O. Theoretical estimate of the orbital pe- 
riod derivative in the low-mass X-ray black hole binary 
XTE J 1 1 18+480 due to magnetic braking and mass loss, 
P mb,ml given in ms/yr, de rived from the second term of eq. 28 
in Johannscn et al. (2009), as a function of two parameters: 
(3 — A / b 1 1 / A /9 : the fraction of mass lost by the secondary 
star that is captured by the black hole; and 7, an index that 
characterizes the strength of the magnetic braking. The spe- 
cific angular momentum, j w , carried away by the mass lost 
from the system has been set to 1 which indicate that all the 
mass lost from the system is lost from the neighborhood of the 
black hole (or its accretion disk, sec IPodsiadlowski et al.120021) . 
Dashed lines represent the value derived from observations and 
the dotted lines, the values at 1-er ver the uncertainty on the 
orbital period derivative, AP. 



star in XTE Jl 118+480 descends from a mag- 
netically peculiar Ap/Bp star that has retained 
most of its primordial magnetic field. In such 
stars very high magnetic fields of ~ 20 kG have 



been observed ( Elkin et all 120101 ) . Binaries with 
intermediate-mass Ap/Bp stars have indeed been 
postulated as the progenitors of compa ct black 
hole X-ray binaries ()Justham et al. 20061) . In par- 
ticular, magnetically coupled, irradiation driven 
stellar winds can lead to s ubstantial loss of sys- 



temic angular momentum (jJustham et al.l [2006) 



required to form low-mass X- ray binaries out from 



intermediate-mass binaries ( Podsiadlowski et al 



2002). The detection of CNO-processed mate- 



rial in XTE Jl 118+480 strongly suggests that 
the donor descends fro m an intermediate-mass 
star ( Haswell et al.ll200"2 ). Our discovery of a large 
period decay in XTE Jl 118+480 is also consis- 
tent with this scenario although we cannot rule 
out the possibility that other mechanism, not yet 
identified, plays a major role in the loss of angular 
momentum and hence the evolution of black hole 
X-ray binaries. In any case, our observations sug- 
gests a faster evolution and shorter lifetimes than 
previously assumed. This would help to reconcile 
the population number of low-mass X-ray binaries 
and millisecond pulsar s, a long-standing problem 
in galactic astronomy ([Podsiadlowski et al.ll2002j) . 
Follow-up spectroscopy in the coming years may 
provide a determination of the second derivative 
of the orbital period with strong implications on 
our knowledge on the formation and evolution of 
black hole X-ray binaries. 

The black hole in XTE J1118+480 was possi- 
bly formed in a violent supernova explosion that 
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launc hed the system via a n asymmetric natal 
kick (Guala ndris et al. 2005) from its formation 
region in the Galactic thin disk (jGonzalez Hernandez 
2006) to its present loc ation in the Galactic 
halo (jMirabel et al.ll200ll) . A lower limit to the 
age of the system of > 11 Myr was derived from 
its peculiar location and kinematics. This limit 
has been used to set constraints on the rate at 
which black holes can evaporate in the Anti-de 
Sitter (AdS) braneworld Randall-Sundrum grav- 
ity model (for detail s see the introdu ction in 
Johannsen et a l. 2009; Johannsenl 2009) via the 
emission of a large number of conformal field the- 
ory (CFT) modes (|Emparan et alJl2003h . 

An upper-limit to the asymptotic AdS curva- 
ture radius in t he extra dim ensions was established 
at L < 80 fim ( Psaltid 2007 ). restricting deviations 
from the gravitational inverse square law to man- 
ifest only at distances smaller than L. 

The orbital period evolution in this system was 
studied before but with a less number of measure- 
ments of T n spread in a shorter interval of time, 
thus providin g a period deriv ative consistent with 
zero at 1-a (jJohannsen 2009), and an upper-limit 
of L < 97 /im. In their eq. 2, these authors as- 
sumed no angular momentum loss due to mass lost 
from the system (j w = 0), no accretion onto the 
black hole (/3 = 0), and the parameter 7, which 
governs the strength of magnetic braking , equal 
to zero (see also Fig. 1 in |johannsenll200"9[ ). If we 
adopt the black hole and secondary masses given 
in Table [1] our determination of the orbital period 
decay provides a much tighter constraint on the 
asymptotic AdS curvature radius of L < 35 fim 
at 2-er. This limit is more restricted than the 
best current table-top experimen t upper-limit of 
L < 44 [im (jKapner et al.ll2007l) . We note here 
that the size of the extra dimensions has been also 
recently constrained from the age of a black hole in 
a extragalactic globular clust er, placing an upper - 
limit as low as L < 3 fim (see Gnedin et al.| [2009) . 
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